Exact Finite Size Study of the 
2dOCP at r = 4 and T = 6| 



G. TellezQ 

Ecole Normale Superieure de Lyon, Laboratoire de Physique, 

Unite Mixte de Recherche 5672 du Centre National de la Recherche Scientifique, 

46, allee d'ltalie, 69364 Lyon Cedex 07, France 

P.J.ForresterQ 

Department of Mathematics and Statistics, University of Melbourne, 
Parkville, Victoria 3052, Australia 

An exact numerical study is undertaken into the finite N calculation of the free 
energy and distribution functions for the two-dimensional one-component plasma. 
Both disk and sphere geometries are considered, with the coupling T set equal to 4 
and 6. Extrapolation of our data for the free energy is consistent with the existence 
of a universal term log A, where x denotes the Euler characteristic of the surface, 
as predicted theoretically. The exact finite A density profile is shown to give poor 
agreement with the contact theorem relating the density at contact and potential 
drop to the pressure in the thermodynamic limit. This is understood theoretically 
via a known finite A version of the contact theorem. Furthermore, the ideas behind 
the derivation of the latter result are extended to give a sum rule for the second 
moment of the pair correlation in the finite disk, which in the thermodynamic limit 
converges to the Stillinger-Lovett result. 



1 Introduction 

The two-dimensional one-component plasma (2dOCP) is a model in classical statistical 
mechanics which consists of A^ mobile point particles of charge q interacting on a surface 
with uniform neutralizing background charge density. The pair potential $(r, r ') between 
particles is the solution of the Poisson equation on the particular surface. In the plane 

$(f,r") = -log (|f-r"|/0, (1-1) 



where / is some arbitrary length scale which will henceforth be set to unity. With the 
potential ( p..l| ) and a uniform background of charge density —pb inside a disk of radius R 
{pb = N/ttR"^) the corresponding Boltzmann factor, which consists of the particle-particle, 
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particle-background and background-background interaction, is given by 

where F := q^/ksT is the couphng. We remark that with r/2 an odd integer, (|1.2| ) is 
proportional to the absolute value squared of the celebrated Laughlin trial wave function 
for the fractional quantum Hall effect |[Lau83|| . 



At the analytic level our knowledge of the properties of the 2dOCP comes from two 
main sources. First, for the special coupling F = 2, the exact free energy and correlation 
functions can be calculated for a number of different geometries ||AJ81| , |Cho81| , |Cai81| , 



JT98|| . Second, the 2dOCP is an example of a Coulomb system in its conductive phase 



and as such should obey a number of sum rules (see e.g. |[Mar88||) which typically represent 



universal properties of such a system. We remark also that the exact solutions at F = 2 
have been an important source of inspiration to identify universal properties. 

In this paper we develop exact numerical solutions at the special couplings F = 4 and 
F = 6 for values of up to 11 and 9 respectively. By undertaking this study we are able 
to test the prediction of Jancovici et al. [|JMP94|] that the expression for the free energy 



F as a function of the number of particles A^ be of the form 

(3F = AN + BN^''^ + ^\ogN + ■ ■ ■ , (1.3) 

where x denotes the Euler characteristic of the surface (x = 1 for a disk, x = 2 for a 
sphere). Furthermore we are able to investigate the rate of convergence of the one and two 
point correlation to their thermodynamic values, as well as the accuracy of certain sum 
rules in the finite system. In fact the latter line of investigation leads us to a new sum rule 
valid for general v dimensional multicomponent Coulomb systems in a spherical domain, 
which relates to the second moment of the density- charge correlation function in the finite 
system. We recall (see e.g. | |[Vlar88| ]) that in the infinite system the second moment of the 



charge-charge correlation function is of a universal form known as the Stillinger-Lovett 
condition. Indeed our sum rule ( [4.24[ ) below gives the finite size correction to this universal 



form in systems with a background. 

As an outline of the paper, we note here that in Section 2 formulas are presented 
specifying the partition function and one and two point distribution functions for the disk 
and sphere geometries, with the coupling an even integer, in terms of certain expansion 
coefficients. These expansion coefficients are in general computationally expensive, but 
reasonably efficient algorithms exist in the literature applicable to the cases F = 4 and 
6. Our numerical results our presented in Section 3. The new sum rules are derived and 
discussed in Section 4, while Section 5 concludes with a summary. 

2 Formalism 

Our interest is in the exact numerical computation of the partition function and one and 
two-point correlation functions for the 2dOCP in a disk and on the surface of a sphere. 
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In the former system the Boltzmann factor is given by Two versions of this model 

will be considered: one in which the particles are confined to a disk of radius R (the same 
disk which contains the smeared out neutralizing background), and the other in which the 
particles are can move throughout the plane. These will be referred to as the hard disk and 
soft disk respectively. In the latter system the Boltzmann factor ( |1.1|) is assumed valid also 
for |rj| > -R, even though the one body potential 7rpfe|rjp/2 is not the correct potential for 
the coupling between a particle and the background in this region (according to Newton's 
theorem outside the disk the background creates the same potential as a charge — at 
the origin, so the correct Coulomb potential outside the disk is A^log |rj|). 

On the surface of the sphere the Boltzmann factor is given by 



1 \ 

) e^^^/^ n \u,v, - u,v,\^ , (2.1) 

^<j<k<N 



2Rj 



where u := cos(^/2)e*'^''^, v := —ism{9/2)e~'^^^'^ are the Cayley-Klein parameters and 
(6*, 0) are the usual spherical coordinates. For our purpose it is convenient to consider 
the stereographic projection of this system from the south pole of the sphere to the plane 
tangent to the north pole. This is specified by the equation 



z = 2Re''^tan-, z = x + iy. (2.2) 
2 

We then have 

2RJ ' ' \2RJ 



l<j<k<N 

^ 1 Z--Zk'^ 



2R 



2.1 The cases T = 4p 

For r = 4p, integrals over the Boltzmann factors ( p. . 1|) and ( p.3|) can be performed from 
knowledge of the coefficients in the expansion 



n {zu-z,fP = Y.ci^\2p)m,{z,,...,z^) (2.4) 



where = (/xi, . . . , iin) is a partition of pN{N — 1) such that 

2p{N -1) > I2i > ■ ■ ■ I2N > 



and 



is the corresponding monomial symmetric function (the rrii denote the frequency of the 
integer i in the partition). The key point for the utility of ( |2.4| ) is that with zj = rje^^^ 
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the are orthogonal with respect to angular integrations: 

roo p27T r27v 

dririg{rl)--- drNTNgir^) d9i--- dO^ m^{zi, . . . , ZN)m^{zi, . . . , zn) 
Jo JO Jo 

= ^,^Tr—yi[G,, (2.5) 

where G^, := 2 Jl^ drr^^'^^'g{r'^) for arbitrary g{r'^). Thus, after also noting that 

n kfc - ^.f' = IK^fc - ^^f'Xli'^k - z,f'. (2.6) 

j<k j<k j<k 

we see that for T = 4p 

J^^ JR^ j<fc ^ lUmil ^^-^ 

(2.7) 

In the case p = 1 this formahsm has been utihzed by Samaj et ah |5PK94], who fur- 



thermore presented an algorithm for the computation of {c^} in this case. Let us now 
consider this latter point. 

In general the coefficients c^^-* (2p) can be calculated from the formula 



c(f)(2p) = — ^ d9i e-'"'^' ■■■ r c/^^ve"^^^^^ \{ (6*^*= - e'^^f^, (2.8 

(ZTTj JO Jo j^j^ 



which follows from (|2.4| ). Since we require \n\ = pN(N — 1), the integral over 9n can be 
performed by changing variables 9j 9j + Ojy {j = 1, . . . , N — 1) to give 



2n r'in 



c(f )(2p) = [ " dOi e-'^^'^ ■■■ r dON-ie-'^'^-''^-' I] - ^"'"f" 

yZ:lX J Jo J j 

X n (e"'^" -e'^^yp. (2.9) 

l<j<k<N-~l 



The simplest case is = 2, when the sum over pairs in (|2.9| ) is not present. Expanding 
(1 — e*^^)^^ according to the binomial theorem gives 

where /xi = p, p + 1, . . . , 2p (for /ii = p we have /xi = yU2 and thus m^^ = 2, while in all 
other cases fii ^ fi2 and so m^^ = m^^ = 1). Substituting in ( p.7|) we see, after some 
minor manipulation, that 

/2p^^ 



/ drigirh [ df2g(rl)\r2-ri\^^ = n'^Y^i ^) / drr'^gir) [ drr'^^ ^gij)- 
JR2 J-R? ^ \ 11 J Jo Jo 

(2.10) 



4 



To calculate c^^\2p) via this method for a general value of would require expanding 
^{N — 1)N products via the binomial theorem, giving a total of (|(iV — 1)N)^p^^ terms 
to determine each value of c^. Thus for a given value of the complexity increases 
exponentially with the coupling p. As we want to determine the for a sequence of 
values of A^ as large as possible, we are therefore restricted to the case p = 1. 

In fact the case p = 1 allows ( |2.8| ) to be computed without using the binomial expansion 
SPK94|]. Instead one uses the Vandermonde formula for the product of differences as a 



determinant to expand the products in ( p.8| ). This gives 

N N 
Pe5jv QgSjv k=l PgSn QeSN k=l 

(2.11) 

which is the formula we used to compute our data in the case p = 1 for A^ = 3, . . . , 10. 



2.2 The cases r = 4j3 + 2 

With r = 4p + 2, decomposing the product of differences analogous to ( |2.6|) shows that 
we must consider the product of differences raised to an odd power. The analogue of ( |2.4| ) 
is then the expansion 

n (zk - z.f'^' = E cr(2P + +^-'^2^^+^-' ■■■zT) (2.12) 

where 2p{N - 1) > /xi > /i2 > ■ ■ ■ > /iiv > 0, EjLi /^j = pA^(A^ - 1) and A denotes 
antisymmetrization. Factoring out the antisymmetric factor Y\j^k{zk — Zj) from both 
sides then gives 

n iz,-z,y^ = J2i''\2p+l)S,{z„...,z^) (2.13) 

^<j<k<N At 

where S"^ denotes the Schur polynomial indexed by the partition /x. Furthermore, analo- 
gous to the orthogonality ( p.5|) we have 

r2TT r^TT 



/ ■ ■ ■ / Ildri rig{rf) / dOi ■ ■ ■ d6N TT l^i ~ Sf,{zi, . . . , ZN)Sf,{zi, 
Jo Jo fj, Jo Jo 



zn) 



Thus for r = 4p + 2, instead of ( |2.7| ) we have 



N 

(5^,,Ar!7r^n^M.+A^-i- (2-14) 
1=1 



N 



I^M = Ar!7r^^(cf )(2p+ l))2^G'.,+7v-^ (2.15) 



ij- 1=1 



According to ( |2.12| ) the coefficients c^^\2p + 1) can be computed from the formula 
with fij + N — j and 2p \—>- 2p + 1 , or equivalently (p.9|) with the same replacements 
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In the case N = 2 this latter formula gives 



c(f)(2p+l) 



-1) 



'2p + r 

Ui + l 



with yUi = p, . . . , 2p. This in turn implies that the formula ( p^.lOD again holds with 2p 
2p+l. 

To obtain data for consecutive values of A^, the computationally simplest case is p = 1. 
However algorithms based on ( |2.8|) (with fij i— > fij-\-N — j and 2p h-^ 2p+l) are inferior to 
methods that determine c\f'^{3) from (|2l3| ) ||FGIL94| , |Dun94| , |STW94|| . The most efficient 
algorithm appears to be the one of Scharf et al. |^TW94|| , where the coefficients c|j^^(3) 
are determined up to = 9. Fortunately the authors of ||STW94|| have kindly supplied 



us with their data (up to = 8), so we do not need to repeat the calculation. 



2.3 The sphere 

The Boltzmann factor for the sphere, stereographically projected onto the plane, is given 
by the r.h.s. of ( p.3|) . Thus, with fj '^Rr j we require 

^(r2) = (i + ^2)-(^-i)r/2-2 ^2.16) 

in the integral (|2.7|) . However, computational savings can be obtained by first noting that 
because the sphere is homogeneous, one particle can be fixed at the north pole, reducing 
the number of integrals from to — 1 (we must also multiply by vr - the area of the 
surface of a sphere of radius 1/2). Thus we have 



N ^ 

w^'^"""^'^n n + b.|2)(7v-i)r/2+2 n i^fc-^^'i 

•"y^ I i=i \^ ^ ) l<j<k<N 



IT 



N~l 





\Zi 


r 


(1 + 




2~)(Af-l)r/2+2 



■ --d^N-i n n , u i2Viv-i)r/2+2 n i^fc - ^.-r, (2.17) 



l<j<k<N~l 



and so should choose 



9^^') = (i + ,2)(^-i)r/2+2 (2-18) 
in (^). 

With gir"^) given by ( p^.lS] ), the formulas ( PT^ ) and ( |2.15| ) show that at F = 4 and 
F = 6 the canonical partition function 

1 



'NT 



dfi - ■ ■ dfN e 



A^! J(R2)'' 

can be represented by the series 

,(JV+1)2 V+l -I N 



^^+1,4 - iv+1 ^V^'^ ^^^) YlW.fA (2Ar+l)! ^^-^^^ 
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V 

To obtain these formulas use has been made of the definite integral 

■dr= ^ '-. (2.21) 



oo 



(l + r)9 r(g) 

Because the sphere is homogeneous, the two-point distribution p(2)((6',0), {9',(j)')) can 
be computed with one particle at the north pole {6' = say). We then have 

P(2)((^^,0), (e',0')) = P(2)W 

SO the two-point function can be computed from an integral of the form ( p.7| ). In fact 
with g{r'^) given by ( 2.16| ) we have 



P(2)(^) = 77^7 — TT 1™ — — ; — (1 + a; ) (1 + a; ) ' \ (2.22) 

where x = tan(^/2). For F = 4 this gives 
o (2iV-l)! 



P(2)(^^)-P,^2(l+^2)2;V-2 



1 ^2/ife 

E (4''H2))^T7TrTn/^^K2iV-2-/i„)!E 



E,(cr(2))2nf^ /i.!(2iV - 2 - p,)! 
while for F = 6 we deduce that 



('3j\r _ 2V 

m(g)=P^ ^2(l+^2)3;V-3 (3^-2) 

E (c(f)(3))^n(^^+^-0!(2iV-3-/i.+.)!E 



^ .,..=0 (/i.+iV-fc)!(2iV-3-p,+A:)! 

E^(i^^(3))2 nili(p. + iV - z)!(2iV - 3 - /i, + ^)! 

2.4 The disk 

In the case of the disk, (|1.2| ) with Vj Rfj shows we require 

9{r') = X(r)e-r^l^^^l^/^ (2.25) 

where x = 1 for < 1 and zero otherwise in the case of the hard disk, while x = 1 for 
all r in the case of the soft disk. Thus from ( |2.7| ) we have at F = 4 

^softdis. ^ e^^v2(^_L) vr^x:(cr(2))^(n-^'j n/^^' (2-26) 

^haMdisk ^ e^^^/^(-L) vr^E(cr(2))^T7^n7(p. + l,2iV) (2.27) 
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while at r = 6 use of ( |2.15|) gives 

N 

Z^f = p^/2jV-3^V23-^(37V-l)/2^37V/2g9^V4^(c!f)(3))2 J] l{l^k + N - k+ I, 3N) (2.28 



M k=l 



with the soft disk case obtained by replacing the incomplete gamma functions by complete 
gamma functions. 

Unlike the situation with the sphere, the density is a non-constant function in the disk 
geometry. Now, with 5'(r^) given by ( |2.18|) we have 



At r = 4 this gives 

P(i)(r) = 2p,e-2-^^^' 



1 ^ ^ (97r Di.r'^)!^!' 

E(^r(2))^77TrTn7(/^. + l,2iV)E ^ 



^,(<^f\'^)?Y^. nf=i 7(/^, + 1, 2iV) 

while at r = 6 one obtains 



El^'^Hs))^ n 7(/i, + AT - , + 1, 3iv) ^ -i:^ 



x-^^ ' " : 1 (2 30) 

E,(cr(3))2 nf=i 7(/i, + iV - J + 1, SAT) 

The corresponding formulas for the soft disk are obtained by replacing the incomplete 
gamma functions by complete gamma functions. 

Finally, we consider the two-point function in the disk geometry. In general this 
quantity is not just a function of the distance between particles, and so we cannot use 
the formalism based on the orthogonalities ( p.5| ) and ( |2.14|) . However, with one of the 
particles fixed at the origin (r*' = say) we have P(2){'^,'f^') = P(2)(^); so in this case 
the formalism used to compute the densities can again be used. Thus using the general 
formula 

we find for the hard disk case 



; 7(pfc + N-k + l,3N) 



P(2)(r) = 4p^e-2w^ 

(2.31) 



E,,,..o(c(^n2))^j^ m <N + 1, 2Ar) EtY ^gg&v) 



SM(ci''H2))2^nf=i7(/i. + l,2Ar) 
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/^(2)(r) = 9pge-3w 



E,(cr(3)P nf=i 7(/i, + iV - J + 1, 3iV) 

for r = 4 and F = 6 respectively. Again the corresponding results for the soft disk are 
obtained by replacing the incomplete gamma functions by complete gamma functions. 



3 Numerical results 



3.1 Free energy — sphere geometry 



In the Introduction it was commented that the free energy is expected to have a large 
expansion of the form ( |1.3|) with x = 2 in sphere geometry. In fact the constant B in 
( |1.3| ), which is a surface free energy, should be identically zero in sphere geometry, so we 
expect a large A^ expansion of the form 



/3F = AA^+-logA^ + C + -- 
6 



(3.1) 



As noted by Jancovici et al. ||JMP94|| , the validity of ( p.l|) can be explicitly demonstrated 
at /5 = 2 because of an exact solution due to Caillol [CaiSl]. The mechanism for the 



exact solution can be seen within the present formalism. Thus, at F = 2 we require the 
coefficients c^^^(l) in ( |2.12| ). But this follows from the Vandermonde expansion (recall 
( 2.11| )), which gives c^^\l) = 1 for /i = 0^ and c^^\l) = otherwise. Substituting in 
( |2.15| ) with (7(r^) given by ( |2.18| ), and making use of ( |2.21| ) we thus obtain ||Uai81|| 



^sphere 



TT 



k=l 



This substituted into the general formula 



NT 



\ogZ 



NT 



(3.2) 



(3.3) 



leads to the expansion | JMP94 



/3F = AT/ + 1 log AT + 1 - 2C'(-1) + o(l) 

D iZ 



(3.4) 



where / = | log(pb/27r^). We remark that by introducing the Barnes G function according 
to 

G{z + l) = T{z)G{z), G(l) = l 



we can write 



N 



X{{k~l)\=G{N + l). 



k=l 
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m + 1) ~ — log AT - -iV^ + - log 2vr - - log iV + C'(-l) " + O(l^) • (3-5) 



The large N expansion of the Barnes G function is known to be |ParOO | 

r2 

This together with Stirling's formula allows us to extend ( p.4|) to the expansion 

= iV/ + 1 log iV + i - 2C'(-1) + ^ + O(^) . (3.6) 

In the cases F = 4 and F = 6, by following the numerical procedure detailed in the 
previous section, we have been able to compute the partition functions (|2.19| ) and (|2.20|) 
up to 11 and 9 particles respectively. The results are listed in Table 0. Our results are 
presented in decimal form. However the terms in the summations of ( [^.19D and ( |2.20D 



are all rational numbers, and we have also calculated the sum itself as a rational number. 
A point of interest is the factorization of the denominator and numerator of the rational 
number. The exact result ( p.2|) shows that at F = 2 only small integers occur in this 
factorization. However our exact data shows that this feature is no longer true at F = 4 



or F = 6. For example, at F = 4 and with = 9 we find that the summation in ( |2.19| ) is 
given by the ratio of primes 

19 • 23 ■ 31 ■ 404431651134013 ■ 56827 
2283125778118138178 • 

To analyze our data we first sought fitting sets of consecutive values of to the ansatz 

/?Fr,7v = ^rA^ + Kr log + Cv (3.7) 

The results are contained in Table |^. Notice that at F = 4 the value of the free energy 
per particle A appears to have converged to 3 decimal place accuracy, while the value of 
K appears similarly to be converging, with the final value in the table differing from 1/6 
only in the third decimal. The general trends are the same for the F = 6 data, although 
the convergence rate (as determined by the difference between sequential values) is slower. 

Next we sought fitting four consecutive values of A^ to the ansatz 

pFr,N = ArN + Kr log N + Cr + Dr/N. (3.8) 

The results of this fit are presented in Table ^. At F = 4 this markedly improves the 
convergence rate, with the final estimate of K now differing from 1/6 by only 3 parts 
in 10^. However at F = 6 the convergence rate is in fact worsened, indicating some 
illconditioning when the extra free parameter is introduced. Note also that the coefficient 
of 1/A^ in both cases appears to be non-zero, as distinct from the situation at F = 2 
exhibited by the analytic result (|3.6|) . 

Finally, we sought to estimate from our data an accurate as possible value of the free 
energy per particle, /3/r say. For this purpose we fitted the data to the ansatz 

1 f ErlN^ F = 4 

l3Fr,N = ArN+-\ogN + Cr + DrlN+\ J' ' (3.9) 
D 0, 1=6 
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Table 1: Exact numerical computation of the expressions ( |2.19|) and ( p.20| ) (in the latter 
case we have set pb = 1), and the corresponding free energy (|3.3|). 







/5i^iV,4 


3 


9.770695753081390794542103296367E+02 


-6.884557862719257767291929292830 


4 


1.081868103379375397165672403770E+04 


-9.289029644211538110263324038604 


5 


1.209528877878741526102013133936E+05 


-11.70315639163470461293716934684 


6 


1.360835037494310939624360869217E+06 


-14.12360906745006986750189927991 


7 


1.537846289459171693753614603094E+07 


-16.54847857521316551691816164401 


8 


1.743564157878398325393942744018E+08 


-18.97661212873318180330363390257 


9 


1.981770773388678655915061613417E+09 


-21.40725661197234419004446417460 


10 


2.257011016434890100740949944465E+10 


-23.83989230877186989649422160272 


11 


2.574639922522006241714385546434E+11 


-26.27414571135846506646694529338 




Zn,q 




2 


781.80154948970530457541038293910180 


-6.661600935308419284761353568226471 


3 


24731.016946702464115291740435512837 


-10.115813481655518642906626162676076 


4 


798906.45662411908447403801186279894 


-13.590999142330226359670889161470696 


5 


25990836.664099377843271224794515169 


-17.073254597869416657276355484106596 


6 


851167572.30792422833993160492670601 


-20.562119579383207945093207167461793 


7 


27989023411.960800446597844273994987 


-24.055078249259894430456119939885817 


8 


923260788226.64381072982338145761830 


-27.551177575665397081224942401207047 


9 


30529687045074.352434196537904510620 


-31.049720671888250916196597607309575 



Table 2: Fitting the values of (5Fy^n with N as specified, taken from Table ^ to the 
ansatz ( ^l7|) . 



N 


^4 




C4 


A, 




c. 


3,4,5 


-2.447509 


0.149600 


0.293616 


-3.526411 


0.178065 


0.267797 


4,5,6 


-2.448705 


0.154963 


0.290968 


-3.506699 


0.109543 


0.283938 


5,6,7 


-2.449038 


0.156787 


0.289696 


-3.515359 


0.145316 


0.269664 


6,7,8 


-2.449271 


0.158300 


0.288384 


-3.516438 


0.152316 


0.263596 


7,8,9 


-2.449423 


0.159440 


0.287231 


-3.516820 


0.155176 


0.260704 


8,9,10 


-2.449524 


0.160290 


0.286264 








9,10,11 


-2.449594 


0.160960 


0.285428 
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Table 3: Fitting the values of PFr^N with N as specified, taken from Table ^ to the 
ansatz (|3.8|). 



N 




Ki 




L»4 


Ae 


Ke 


Cg 


De 


3,4,5,6 


-2.450743 


0.175200 


0.258672 


0.049566 


-3.5382 


0.3594 


0.0572 


0.4839 


4,5,6,7 


-2.449773 


0.165568 


0.2740449 


0.025973 


-3.5086 


0.0654 


0.4121 


0.2363 


5,6,7,8 


-2.449905 


0.167146 


0.2712323 


0.031065 


-3.5193 


0.1932 


0.1842 


0.1417 


6,7,8,9 


-2.449914 


0.167268 


0.2709949 


0.031065 


-3.5180 


0.1748 


0.2199 


0.0779 


7 ,8,9,10 


-2.449896 


0.166989 


0.2715743 


0.029956 










8,9,10,11 


-2.449892 


0.166917 


0.2717321 


0.029634 











Table 4: Fitting the values of (3Fy^n with N as specified, taken from Table ^ to the 
ansatz 



N 


Ai 


Ca 


Di 


Ea 


^6 




1)6 


3,4,5,(6) 


-2.4501031 


0.275576 


0.012460 


0.026276 


-3.513916 


0.205966 


0.110598 


4,5,6,(7) 


-2.4498406 


0.271639 


0.031880 


0.005215 


-3.518863 


0.250494 


0.011648 


5,6,7,(8) 


-2.4498809 


0.272364 


0.027574 


0.003235 


-3.517146 


0.231609 


0.063153 


6,7,8,(9) 


-2.4498875 


0.272503 


0.026605 


0.005465 


-3.517466 


0.235770 


0.049709 


7,8,9,(10) 


-2.4498842 


0.272423 


0.027240 


0.003788 


-3.517540 


0.236870 


0.045600 


8,9,10,11 


-2.4498841 


0.272420 


0.027272 


0.003695 









thus assuming the universal term in ( |3.1| ). Four free parameters are used at F = 4, while 
only 3 free parameter are used at F = 6, in keeping with observed illconditioning when 
a fourth parameter is introduced. Our results are presented in Table ^ where /3/r is 
determined by Ay. We see that there at F = 4 we appear to have convergence to 7 digits 
with the estimate 

/3/4 = -2.449884 ■• • (3.10) 
while at F = 6 our final estimate is 



^/g = -3.5175... (3.11) 

accurate to 5 digits. 

We note that there is some early literature on estimating (Sf^ and fif^ from exact small 
N numerical data ||JM83|| . Using only the values of PF^^r for = 1, 2 and 3, the quantity 



3F F 

^/r = ^/r + ( Y + l) + 4 log vrp6 - log 

was estimated for F = 4, 6, . . . , 10. In particular, at F = 4 and F = 6 these estimates of 
Pfr give 

/3/4 ^ -2.1585, pfe ^ -3.330, 
which differ from our estimates (|3.10|) and ( p.ll|) in the first decimal place. 
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3.2 Free energy — disk geometry 



For the disk geometry, the prediction (O) gives a large N expansion of the form 

f3Fr = AN + BN^/^ + — logN + C + --- (3.12) 



As in the case of the sphere geometry, this prediction can be verified analytically using 
the exact solution for the isotherm F = 2 | AJ81| |. The exact solution gives [JiVlP94 

^^hard _ ^^^^ ^ /?72iVl/2 + _L log jV + 0(1) (3.13) 

where 

f3f2 = -\og{p,/27T'), Pl2 = V2j^ dy log {-{1 + erf y)). 

Some details of the expansion of (3F2 are different for the soft edge version of the OOP in 
a disk (recall Section 1). From the exact formula 



7Soft 
'N,2 



IT 



and the asymptotic expansion ( p.5|) we see that 



soft 



/3/iV + -logiV-C'(-l) 



1 



O 



(3.14) 



720A^2 

Thus indeed both ( p.l|) and ( 3.14 ) contain the universal term (1/12) log A^, although (|3.14 ) 
does not contain a surface tension term (this fact has been noted previously in [[FGIL94|| ). 

At F = 4 and F = 6 we obtained exact numerical evaluation of the partition functions 
(12:2^ ), (|]23) and (|2]28|) (and the modification of (|]28D for the soft disk case) as in the 
sphere case. Our results for the corresponding value of (3F are contained in Table ^. To 
test the prediction (p.l|), we sought to fit our data to the ansatz 



PFn,v = /3frN + BrN^'^ + log A^ + Cr + | 



soft disk 
hard disk 



(3.15) 



where /?/r is given by ( |3.10| ) and ( p.ll| ) for F = 4 and F = 6 respectively, and the choice 
in ( |3.15| ) is made retrospectively on the criterium of obtaining better convergence. 

Our results are obtained in Table ^. We see that for the hard disk at F = 4 our final 
estimate of A'4 differs from 1/12 by only 2 x 10^'*. At F = 6 we see that more data would 
be needed to get a stable sequence, although the final estimates of Kq are consistent with 
the expected value of 1/12. 



3.3 Density and two-point distribution 
Density 

Consider for definiteness the disk geometry with a hard wall at F = 4. Using the formula 
( |2.29| ) the density profile can be calculated for up to 10 particles. One way to present the 
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Table 5: Exact decimal expansion of the free energy for the hard and soft disk at F = 4 



and r = 6 




N 




psoft 


3 
4 
5 
6 
7 
8 

y 

10 


-6.07705853011644579848828232852953 
-8.30894530308837749094468707356467 
-10.5685824419856069054748395707000 
-12.8480499008173510151377678768908 
-15.1423987396644292302500680775824 
-17.4483520149155330139161065965798 
-19. mo6864y04052iz 10598150968 /4zio 
-22.0868149972503557220763154840028 


-6.38430353764202167882687100789504 
-8.6/246/ / 1929039/ 19253598118439664 
-10.9817913623032469741300225072724 
-13.3060582270200975291371029447052 
-15.6414978836634761215222474874096 
-17.9856458201720068377211714643235 
-20.3368227969313363262711724430690 
99 fiQ'^897SQ7'ifi97nn'^'^'^fi9R'^88n871 ^A'\ 


N 






3 
4 
5 
6 
7 
8 


-9.0582041809587470427592556776938317 
-12.306265058620940233015626198772823 
-15.583591141405785588643527765993475 
-18.886678348734296934648840469575921 
-22.209056127812161704085770192533417 
-25.545482070626355796809539664033139 


-9.1916690110088058948684895153913657 
-12.467150515773535356614120708869630 
-15.769625685129047660199805300971936 
-19.095091912250933709000748332754870 
-22.437790137971372572352358156488860 
-25.793196864919170855940747024418097 



Table 6: Fit of the ansatz ( p.l5| ) to the data of Table 



N 


^hard 


^hard 


^hard 


^soft 


^hard 


^soft 


^soft 


3,4,5,(6) 


0.749371 


0.059801 


-0.091054 


0.497409 


0.120202 


-0.052807 


0.073687 


4,5,6,(7) 


0.728988 


0.081365 


-0.080181 


0.509625 


0.099801 


-0.040616 


0.040317 


5,6,7,(8) 


0.723951 


0.087261 


-0.078408 


0.522124 


0.076905 


-0.022675 


-0.004884 


6,7,8, (9) 


0.726340 


0.084219 


-0.078810 


0.521397 


0.078345 


-0.024029 


-0.001556 


7,8,9,(10) 


0.727263 


0.082957 


-0.078795 


0.518587 


0.084298 


-0.030427 


0.014202 


8,9,10 


0.726874 


0.083523 


-0.078873 











N 


rjhard 


j^hard 


/^hard 


rjsoft 


TV'SOft 


^soft 


3,4,5 


1.104506 


-0.092158 


-0.317518 


0.967066 


-0.059461 


-0.248851 


4,5,6 


0.884919 


0.140146 


-0.200388 


0.795791 


0.121733 


-0.157491 


5,6,7 


0.874984 


0.151776 


-0.196890 


0.786513 


0.132593 


-0.154224 


6,7,8 


0.951461 


0.054407 


-0.209757 


0.842635 


0.061139 


-0.163667 
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Hard disk density 

r=4 




Figure 1: Density profile in the hard disk case for several values oi N ai T — A. The 
boundary of the disk is taken as origin. 
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data is in graphical form with the boundary of the disk taken as the origin. This is done 
in Figure [^. The plot shows rapid convergence of the profiles near the boundary. 

To investigate the rate of convergence of the whole profile as measured from the 
boundary to the thermodynamic value we can investigate the contact theorem ||CFG80 



This expresses the thermodynamic pressure in terms of the density at contact with the 
wall, and the potential drop across the interface (which in turn is proportional to the first 
moment of the density profile). Explicitly the contact theorem states 

[1 - -jpb = p(0) - 2'KpbT x{p{x) - pb) dx (3.16) 

where we stress again that the density is measured from the boundary. 

Much to our initial surprise, the convergence of the r.h.s. to the l.h.s. for the finite N 
data is very slow. For 10 particles the error is of order 30%. Further investigation reveals 
that this is not special to the coupling F = 4. At F = 2 we have the analytic expression 



JanSl 



p^.r = — ; , 0<r<R 

where r is measured from the boundary and the background density is taken to equal 
unity. Choosing = 10 and substituting in (|3.16|) again gives an error of order 30%. 
Indeed choosing A^ = 500 still gives an error of order 3%. 

In fact the slow convergence of ( p. 16 ) can be understood analytically by making use 
of a sum rule for the OCP applicable for the finite disk |PFG80|| . This sum rule reads 

P(i)(0) - (l - \)p, = C'-'ipW^^ - ^) - P^) (3.17) 

where P(i){r) is measured inward from the boundary. Noting that charge neutrality re- 
quires 

r (p(i) (R-r) - pb)dr = {R - r) (p(i) (r) - pf.) dr 



we can write 



N 



I r\p{i){R - r) - Pb 
TpW 



N 



dr 

{R-r)\p{i){r)- Pb) dr 
£\-2rR^ + ?>r'^R - r^)(p(i)(r) - pb) dr 



N 

/ \ 3F(p;,7r)'^/2 f \ 

SrpfeTT r(^p(i)(r) - Pb) dr — ^%P(i)(^) - Pb) dr 

^^(P(I)W - Pb)dr 



N 



This shows that the finite size corrections to the r.h.s. of ( |3.17| ) are proportional to A^ ^Z^, 
thus explaining our empirical observation. 
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Sphere two-point correlations 

N=8 




Figure 2: Two-point correlation in the sphere case for = 8 particles at F = 4 and 
F = 6. 

Two-point function 

At F = 2 and in the thermodynamic limit the two-particle distribution function has the 
exact evaluation ||J an81|| 

This is a monotonic function, with the corresponding truncated distribution p^)(0,r) : = 
P(2)(0,'r) — P(i)(0)p(i)(r ) exhibiting Gaussian decay to zero. There is evidence, both 
analytic and numerical ||J anSl] , |CLWH82| which suggests that for F > 2 the two-particle 



distribution exhibits oscillations. At F = 4 this feature has already been observed in the 
exact finite A^ calculation of p(2)(0,r) by Samaj et al. [^PK94|] . Furthermore, this feature 



should become more pronounced as F increases. This is indeed what we observe when 
plotting our results for F = 4 and F = 6 on the same graph (see Figure 

The fact that the 2dOCP is a Coulomb system in its conductive phase implies that 
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in the bulk the second moment of the truncated distribution obeys the Stilhnger-Lovett 
sum rule 

L^"'V5)(0,r)cir = --. (3.18) 
For the hard disk in the finite system we can compute 



^ P 2)(0,r)c?f (3.19) 

r\<R 

and compare it with the universal value given by ( p.l8| ). At F = 4 and with = 9 we 
find agreement with the universal value to within 2%. In fact, analogous to the integral 
in (|3.17|) , the integral ( p.l9|) can be evaluated exactly and the terms which differ from 
—2 / ttF read off. In the hard wall case we find 



while in the soft wall case the same expression results except that the boundary term 



NpLAO, R)/ pI is no longer present on the r.h.s., while on the l.h.s. the integral is over 



R2. 



We see from ( p.20|) that the deviation in the finite system from the bulk value (|3.18 
is determined by 

-^((P(i)(0) - p,)/p, + Arp5)(0, R)/pI 
and thus consists of a bulk and surface contribution. 



4 New sum rules 



In this section we present the derivation of the sum rule ( p.20| ) and its generalization to 



multicomponent Coulomb systems. First we show that the sum rule can be derived within 
the formalism of Section ^ then we present a more general derivation of the sum rule. 



4.1 The case F even 

The formalism presented in Section |^ is valid only if F is an even integer. Within this 
formalism we can use the expressions ( p.31| ) and ( |2.32| ) for the two-point correlation 
functions (and its generalizations to higher F) to compute the second moment 

J^f^pl^{0,f)df, (4.1) 

where A is a disk of radius R (hard disk) or (soft disk). For example in the hard disk 
case with F = 4, for each term in the sum ( |2.31| ) the integral ( |4.1| ) gives an incomplete 
gamma function 'j{pj + 2, 2A^). Then we use the recurrence relation 

7(/x, + 2, 2N) = (/i, + l)7(/i, + 1, 2N) - e-^"" {2NY^^^ (4.2) 
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to split the expression in two. The first term is proportional to P{i)(0) while the second 
is proportional to p(2)(0,i?). The sum rule (|3.20|) follows from that. 



The calculation can be easily generalized to any even F. In the soft disk case since the 
incomplete gamma functions are replaced by complete gamma functions the recurrence 
relation ( [4.2|) does not have a second term on the r.h.s., therefore there is no surface 
contribution proportional to pj2-j{0,R) in the sum rule. 

4.2 General case 

In fact a more general derivation of this sum rule, valid for any value of the coupling 
constant, can be obtained by studying the variations of the density as a function of the 
size of the disk. 

Let us consider the general case of a multicomponent jellium in u dimensions confined 
in a spherical domain A of radius R and volume V = Vt^R'^ /u with Vt^ = 27r'^/^/r(z//2). 
The system is composed of s different species with charges {ea)ae{i,...,s} and there are Na 
particles of the species a. Let N = J2a be the total number of particles and let us 
define the average density of the species a, pa = Na/V and the total average density 
p = N/V. As in the preceding sections pb is the background number density and let Cf, be 
its charge so that the background charge density is etpb- For convenience let us define the 
"number of particles of the background" by Nf, = p^V . In general the Coulomb potential 
is 

{— In r, if z/ = 2 
, otherwise, (^•"^) 
Z/-2' 

and the Coulomb force is 

F{r) = -V<l>(f ) = ^ . (4.4) 
The Hamiltonian of the Coulomb system is 

1 ^ r e^p^ r 

U = oIIe„^e„.'^(^i-^i)+ebPbIIe„, / df^fi-f) + ^ / dfdf'^f-f'). (4.5) 

2 .^^ Ja 2 Ja2 

We shall consider the correlation functions in the canonical ensemble 

pS...aAru . . . , r .) = ( E ■ • ■ E ^(^"1 - ^""..J ■ • ■ ^i^n - r „„,J ) . (4.6) 

\il=l in = l I 

The (■ ■ ■) is the average in the canonical ensemble and in the preceding sums if some 
aa = Oih we exclude the term ia = % as usual. 

In three dimensions in order to have a well defined thermodynamic limit we shall 
restrict ourselves to the case where all electric charges Cq, have the same sign and the 
background carries a opposite neutralizing charge. In two dimensions we can also consider 
systems with charges of different signs and eventually without background {pb = 0) if the 
coupling contants |/3eQ,e^| < 2 for all pair of charges (cq,, e^) of different signs. 
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4.2.1 Contact theorem sum rule 



The derivation of the sum rule for the second moment of the two-point correlation function 
is similar to that of the contact theorem for a spherical domain [|CFG80|| . Let us first 
show here the generalization of this contact theorem for the multicomponent jellium. We 
consider the canonical partition function (times N\) 



Q 



df'^ exp{—j3U) 



(4.7) 



as a function of the volume V. We shall compute the thermodynamical pressure p*-^-* = 
c^log Q/dV in two different ways. The derivative is done at fixed number of particles and 
fixed Nh. In general using the scaling r = V^^'^f we have 



oV 



dV Q 
where A is a sphere of volume 1. 

A first way to compute the derivative of U is by using the general formula 

^ = f ■ Fir) . 

dv vv ^ ' 



(4.9) 



This gives, together with the definition (^^) of f/, 

/ dfdf'f ■ F{f - f') ^ Caea'P^al' 
+ebPb / dfdf'{f - f') ■ F{f - f') V CaP^^ 



W(f') 



H — elpl f dfdf'{f — f') ■ F{f — f') 
2 Jh? 



(4.10) 



We can transform the preceding expression by using the first equation of the BGY hier- 
archy 

ksTV pa{f) = CaPbCb df F{f -f')p^^\f)+ I df'y^eaea'F{f -f')p^^l,{f,f) . (4.11) 

Ja Ja ^, 

The r.h.s of ( |4.11 ) appears in the first and second lines of ( 4.10| ). Replacing it by the l.h.s 
of (|;TT]) we find 



P + 



uV 



ksT f dfY^f ■Vp^^\f) 

J A „ 

-CbPb I dfdf'f' ■ F{f - f') V eap[ 
Ja2 ^ 

~^nPb^b / dfdf' (r — f') ■ F[f — r ' 

2 JA2 



(4.12) 



The first term of the r.h.s of the preceding equation can be computed by integration by 
parts while the others can be computed using the definition (^.4|) of the Coulomb force 
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F and Newton's theorem. This yields the following expression for the thermodynamical 
pressure 

-^-^j/rr^Y.^J^\r). (4.13) 

The other way to compute the thermodynamical pressure is to use the actual scaling 
properties of the Coulomb potential 



_ I if = 2 

' ^{r), otherwise. 



(4.14) 



uV 

Substituting this expression in ( |4.8| ) gives 



where Q = J2a ^a^a + ^bNh is the total charge of the system. 



Equating the two expressions ([4.13|) and ( |4.15| ) of the thermodynamic pressure we find 



the generalization of the contact theorem 



P + 5M%-T.<^Ip,) + "-J'I3{U) (4.16) 



where g(r ) = ^aPa\^) + ^bPb is the local charge density. 



4.2.2 Density-charge correlation second moment sum rule 

Similar calculations lead to the second moment sum rule for the density-charge truncated 
correlation function Z]/3 e/3p[fj"^(0, r ). Here we consider the quantity 

. Nc 

Q„=/ cif^e-^^^5(f,J, (4.17) 
•^^"^ i=i 

as a function of the volume V . Note that the density of the species a at the center of the 
spherical domain is p2^'*(0) = Qa/ Q- Like in the preceding section we want to compute by 
two different ways the quantity Q~^dQa/ dV . Using the same scaling argument as before 
we have 



Q dV V 
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Using eq. ( [4.9|) and the definition ( [4.5|) of tlie Hamiltonian U we find 



1 dQa 
Q dV 



^ - 1 (1 



V 



/ c?f rff' f ■ F{f - r') ^ ef3e^pfp^{Q, f, r') 

pa 

+ / f ■F(f)e„^e^pg(0,r) 

+efePf, / (if (if' ^ e/3p[fj(0, f ) f • F(f - f ') 

-CfePfe / dfdf'f' ■ F{f - f ') e/3p[fj(0, f ) 



1 



+ -elpl I dfdr'{r—f')-F{f — f')p'^\0) 
Using the second BGY equation 

kBTVjfp^^^iO, r) = epebPb df' F{f - f ')pf^^j(0, f ) 



(4.19) 



+e^e,F(f)pg(0,f: 



+ /df'F(f-f')EwS,(0, 
and then integration by parts 

E//"* ■ v.pS(o,f) = z^i^EpS(o,^) - Kiv - i)pL'^(o) 



(4.20) 



(4.21) 



we can arrange expression ( [4.19| ) to find, after computing explicitly the integrals involving 
F using Newton's theorem, 



Q-W = 



+ ^^p«(0) 



Z/ + 2 



(4.22) 



The second way for computing Q ^dQa/dV is by using directly equation ( ^.14] ) into 
equation ( f4.18| ). This gives. 



1 dQa 
Q dV 



^P«(0) + f... 



+^/?(f/p«(o) 



V 



E4^/3 



pi'Ho) 



(4.23) 



where p«(0) = 5(fi,„) is the microscopic density of o^-particles at the center of the 
domain A. 
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Comparing the two expressions ([4.22|) and (|4.23| ) of Q'^dQa/dV gives a sum rule for 
the second moment of the density of a particles-electric charge correlation function. The 
sum rule takes a nice form by considering the truncated correlation function and making 
use of the contact sum rule ( |4.16|) , 



2u 



'-'^./3(f/p«(0))^ (4.24) 



In the case of the two-dimensional OCP (z/ = 2, s = 1 and = — ei) this is exactly 
the sum rule ( p.20|) announced in the preceding section 



fVf2)(0,r)rfr-= -AL (0)/p, + iVp5)(0,i?)/p^) . (pD 



The sum rule (|4.24|) is in fact a series of s sum rules for the density- charge correlation 
function Y.p^pPap (0)"^) ^ach species a. By taking the sum of these sum rules with 
the factors Cq, we find a sum rule for the charge-charge truncated correlation function 
S{^,r) = Ea,/3e„e^p2i^(0,f), 



2z/ 



[dfr's{o,f) = Ee«pL'Ho) + — i?^Ee«pir(o,/2) 

a ^ a,l3 



+ ^-^HUJ:^^P^o^\0)) . (4.25) 



a 



4.2.3 Thermodynamic limit of the sum rules 
Canonical ensemble 



In order to study the relationship between sum rules ( [4.24| ) and ( |4.25| ) and the Stillinger 



Lovett sum rule, we need to know the behavior of the correlation functions as they ap- 
proach the thermodynamic limit. This behavior is different depending on the ensemble 
used. In this section we continue to work in the canonical ensemble. 

In general we shall suppose that in the thermodynamic limit the system is in a fluid 
and conducting phase. In this case the density becomes uniform in the thermodynamic 
limit p^\0) Pa and 

{p^^\0)uf^{PaUf = 0, (4.26) 

because in the canonical ensemble the density does not fluctuate. 

Let us flrst consider the case of a multicomponent Coulomb system without back- 
ground (in two dimensions with small Coulomb couplings). In that case equation ( ^.24| ) 



becomes 

pW(0) + ^i?^EpiT(0,i?) = 0. (4.27) 
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This equation (|4.27|) give us the behavior of the correlation functions as they approach 
the thermodynamic hmit 

EpLf (0,i?)~-pPa/iV. (4.28) 

7 

This is a generahzation of an aheady known result concerning the existence of 1/iV tails for 
the correlation functions of one component fluids with short range forces |[LP61|| . However, 
for a neutral system taking the sum of equations ( [4.28| ) with the coefficients Cq, show that 
the charge-total density correlation does not have tails, 

^'^Ee^pif (0,i?)-0, (4.29) 

It is likely that a similar behavior exists in the general case (pf, 7^ 0, v = 2,3), 
so it would be difficult to derive from ( [4 .241) partial sum rules for the density-charge 
correlations in the thermodynamic limit because with the 1/N tails, one cannot commute 
the thermodynamic limit with the integration over the space. However, one can conjecture 
that although the density- density correlations have tails, in the conductive phase the 
total density- charge correlations do not have these tails as it is in the case when pf, = 0. If 
this is true, and assuming that the convergence of the charge-charge correlation function 
is uniform (in order to commute the thermodynamic limit with the integration over the 
space) , one can recover the Stillinger-Lovett sum rule from the sum rule (|4.25|) for finite 
systems, 

^ / drr^S{0,f) = -1. (4.30) 

2l/ JR" 

The fact that we recover the Stillinger-Lovett sum rule is of course not a proof of our 
conjecture, but at least it show that our conjecture is not in contradiction with well known 
results. 



Grand canonical ensemble 



For systems with short range forces the correlations functions do not have tails in 
the grand canonical ensemble as they approach the thermodynamic limit [P^P61|| . We will 
show that this is also the case for two-dimensional Coulomb systems with small couplings 
when there is no charged background and assuming this is also the case in general for a 
multicomponent jellium we will discuss the thermodynamic limit of the partial sum rules. 

The partial sum rules ( ^.24| ) obtained before are different in the grand canonical en- 
semble. The grand canonical ensemble is parametrized by the background density pb and 
s — 1 fugacities {2^} used to fix s — 1 average densities p^, the remaining density fixed 
by electroneutrality. The grand canonical version of the sum rules ( [4.24| ) can be obtained 
in a straightforward manner by adapting the calculations of the last section. However 
special care should be taken because of the fiuctuation of the average densities in the 
grand canonical ensemble. These fiuctuations add some extra terms to sum rule ( |4.24| ), 



2iy 



ldfr^T.e,P^f{0,n = (p«(0)) + ^i?^Epir(0,i?) 
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+^^.2^E^^7PL^H0)^^ (4.31) 

To proceed with the discussion of the thermodynamic hmit of this sum rule, we need 
to use a relation that will allow us to simplify the terms on the r.h.s. of equation ([4 .311 ) 
in the thermodynamic limit. This relation reads for z/ = 2 or 3, 

This relation is a consequence of the scaling properties of the Coulomb potential. To 
prove it, let us consider the thermodynamic grand canonical pressure 

{z,}, p,) = Jim V-^ In {z^}, p,, V) , (4.33) 

where S is the grand canonical partition function. Using the scaling properties of the 
Coulomb potential we have for z/ = 3 

/?p(/?, {^J, p,) = X^PPiXP, {A-3/2zJ, A-3p,) , (4.34) 

and for z^ = 2, 

{^J, Pb) = mP, {A-2(i-(^^^/^«^4, A- V) , (4.35) 

for any positive number A. Taking the derivative of these relations with respect to A, then 
putting A = 1 and using the usual thermodynamic relations yields for u = 3, 

,^(ff),_L,^,,g^(f;),l,_,,g. (4.36) 

and for z^ = 2, 

,,.X:(l-f)p. + /5..f . (4.37) 

where (if) is the total internal energy (including the kinetic term). The announced 
relation ( |4.32| ) follows from taking the derivative of (|4.36| ) and ([4.37|) with respect to the 
fugacities. 

As before let us consider first the case pb = (in two dimensions for systems with 
small couplings). Then equation ( 4.31 ) together with equation (|4.32|) shows that the 
grand canonical total density-partial density correlation function does not exhibit any 
1/A^ tails, 

^''EpLT(0'^)^0- (4-38) 

Now if we suppose that in the general case {pb 7^ 0, z/ = 2, 3) this property still holds and 
that the density-charge correlation functions converge uniformly we recover the partial 
sum rules 

^t*>^Ee,.!.f(0.O^§^. (4.39) 
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that have been previously derived by Suttorp and van Wonderen ||SvW87|| in the three di- 
mensional case. These equations also hold for two-dimensional systems. One can recover 
the Stillinger-Lovett sum rule ( [4.30|) by taking the sum of these equations ( |4.39|) with the 
factors Ca and using electroneutrality. Notice that the condition ( [4.38|) on the thermody- 
namic limit of the two-point correlation function when one of the points is in the boundary 
is different from the usual condition needed to prove the Stillinger-Lovett ||MG83|| that 
the correlation function of the infinite system should decay faster than l/r"^"^. 

Notwithstanding the relation of the sum rules ( |4.24| ) and ( ^.31| ) with the Stillinger- 
Lovett sum rule ( [4.3UD , let us stress that for finite systems these sum rules are not screening 
sum rules like the Stillinger-Lovett sum rule since for finite systems the screening of exter- 
nal charges does not exists (because since the total electric charge is conserved, the excess 
of charge can not leak out to infinity like it does in infinite systems). From the derivation 
presented in the previous section it is clear that the new sum rules should be seen more as 
a second order contact theorem rather than a screening sum rule. Futhermore when there 
is no background {pb = 0) the relation with Stillinger-Lovett sum rule disappears because 
the term containing the second moment of the density-charge correlation vanishes. 



5 Summary and conclusion 

Expanding the power of the Vandermonde determinant that appears in the Boltzmann 
factor of the 2dOCP in terms of simple orthogonal polynomials we have been able to 
develop exact numerical solutions for values of the coupling constant F = 4 and F = 6 
for finite systems up to 11 and 9 particles respectively for different kinds of geometry 
(sphere, soft and hard wall disk). With these solutions we have been able to test the 
prediction [|JMP94|] of universal logarithmic finite size corrections to the free energy ( |1.3D . 



Studying the correlation functions has lead us to find a new sum rule ( p. 20 ) similar to the 



Stillinger-Lovett sum rule for finite systems. This sum rule can be derived within the for- 
malism of section |^, but can also be generalized to higher dimension and multicomponent 
jellium systems (eq. ( [4.24| )). 

Further applications of the formalism presented here are the study of surface corre- 
lations which are expected to have a universal behavior at large distances ||Jan95|| . Also 



the formal expressions of the correlations functions ( p.31| ) and (|2.32| ) could eventually be 
used to find higher order sum rules or other general properties. 
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